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Abstract 

Volume estimates of metric balls in manifolds find diverse applications in information and coding 
theory. In this paper, some new results for the volume of a metric ball in unitary group are derived 
via various tools from random matrix theory. The first result is an integral representation of the exact 
volume, which involves a Toeplitz determinant of Bessel functions. The connection to matrix-variate 
hypergeometric functions and Szego’s strong limit theorem lead independently from the finite size formula 
to an asymptotic one. The convergence of the limiting formula is exceptionally fast due to an underlying 
mock-Gaussian behavior. The proposed volume estimate enables simple but accurate analytical evaluation 
of coding-theoretic bounds of unitary codes. In particular, the Gilbert-Varshamov lower bound and the 
Hamming upper bound on cardinality as well as the resulting bounds on code rate and minimum distance 
are derived. Moreover, bounds on the scaling law of code rate are found. Lastly, a closed-form bound on 
diversity sum relevant to unitary space-time codes is obtained, which was only computed numerically in 
literature. 


Index Terms 

Coding-theoretic bounds, random matrix theory, unitary group, volume of metric balls. 

L. Wei and J. Corander are with the Department of Mathematics and Statistics, University of Helsinki, Finland (e-mails: {lu.wei, 
jukka.corander}@helsinki.fi). R.-A. Pitaval is with the Department of Mathematics and Systems Analysis, Aalto University, 
Finland (e-mail: renaud-alexandre.pitaval@aalto.fi). O. Tirkkonen is with the Department of Communications and Networking, 
Aalto University, Finland (e-mail: olav.tirkkonen@aalto.fi). 

This work was presented in part at 2015 IEEE International Symposium on Information Theory. 


June 25, 2015 


DRAFT 



2 


I. Introduction 

Determining the volume of metric balls in Riemannian manifold, in particular unitary group, is the key 
to understand several coding and information theoretical quantities. Performance analysis of unitary space- 
time codes [1-3] requires the knowledge of volume in the unitary group [4,5]. For channel quantizations 
in precoded multi-antenna systems, the characterization of rate-distortion tradeoff is directly related to 
volume calculations in the manifold of interest [6,7]. Estimating fundamental coding bounds such as 
Gilbert-Varshamov and Hamming bounds relies on the volume of the corresponding metric ball [8-11]. 

Despite the need to accurately estimate the volume of metric balls in the unitary group, results in this 
direction are rather limited. Volume estimates have been derived in [9-11] when the radius of metric ball 
is small. In this paper, we study the volume of a metric ball, valid for any radius, in the unitary group with 
chordal distance. The same problem was considered in [4] in the study of diversity sum bounds of unitary 
space-time codes, where the authors relied on numerical integrations to evaluate the volume. The starting 
point of the current paper is that such a numerical step may not be necessary. Specifically, we show that 
the exact volume boils down to an integral involving a Toeplitz determinant with Bessel function entries. 
This representation gives rise to an explicit formula in the simplest case of a two-by-two unitary matrix. 
We then present a simple asymptotic volume formula, which is the main technical contribution of the 
paper. Two distinct paths that led to this asymptotically exact formula have been identified: one is based 
on the connection of the Toeplitz determinant to a hypergeometric function of matrix argument, the other 
directly invokes Szego’s strong limit theorem on Toeplitz determinants. These powerful tools of random 
matrix theory, albeit being subjects of intense studies for several decades, have not been fully utilized 
in the coding theory community. Surprisingly, the limiting volume formula is already quite accurate for 
dimension as small as three. The reason behind the rapid convergence is examined, where it is found 
that the asymptotic formula approaches its limit super-exponentially fast as dimension increases due to 
a discovered mock-Gaussian property. 

As an application, we study some basic coding-theoretic questions of unitary codes via the derived 
asymptotic volume formula. In particular, analytical expressions for the Gilbert-Varshamov and the 
Hamming bounds on codeword cardinality as well as the corresponding bounds on minimum distance of 
a code are derived. We also derive formulas for bounds on the code rate and its scaling law, which are key 
to establish existence results of unitary codes. In addition, a closed-form upper bound on diversity sum 
of unitary space-time codes is obtained. These simple-to-compute analytical results capture the behavior 
of the bounds reasonably accurately, which also lead to useful insights and properties of the bounds. 
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The rest of the paper is organized as follows. In Section II we formulate the problem considered in this 
paper, where we first define volume of a mefric hall in unifary group before staling the coding bounds of 
interests. Section III is devoted to the derivation of exact and asymptotic volume formulas. In Section IV 
the derived analytical results are utilized in the study of coding-theoretic bounds of unitary codes. We 
conclude the main findings of fhis paper in Section V. Proofs of some fhe fechnical resulfs are provided 
in fhe Appendices. 


II. Problem Statement 


A. Volume of a Metric Ball 

Consider a mefric ball around fhe identity element I„ in the n-dimensional unitary group U{n) with 
Euclidean distance (chordal distance) r, 


.B(r) = {UG[/(n) I ||U-I„||F<r}, (1) 

where 11 • | |f is the Frobenius norm topological metric. We consider the invariant Haar measure defining 
a uniform disfribufion on U{n). For any measurable sef S C U{n) and any U G U{n), fhe Haar measure 
satisfies //(U5) = Due fo fhe homogeneify of U{n), fhe characferisfics of fhe ball (1) centered 

af any group elemenf, say A, would be fhe same. 

The eigenvalue decomposifion of U lakes fhe form 

U = S“^ES, (2) 


where S G U{n) and fhe diagonal enfries of E are n complex numbers on fhe unit circle. More 
precisely, S takes values in a flag manifold of equivalence classes of U (n) modulo diagonal matrices of 
the form E. The joint density of the angles 6i is given by [12] 


p{9l, ...,9n) 



l<j<k<n 


( 3 ) 


where —tt < 0, < tt, i = 1,..., n and the constant c = n!(27r)'^. The product ni<j<fc<n “ e*®*" is 
the Jacobian of the transform (2) of Haar measure (dU) to the measure (dS) nj=i ^9j. In the random 
matrix theory literature, the joint density (3) is referred to as circular unitary ensemble [13] and c is the 
corresponding partition function. 

The condition on the distance measure ||U — I„||f < r in (1) is equivalent to {9i/2) < r^/4. 

Thus, the (normalized) volume of the metric ball (1) equals the following n-dimensional integral [4] 


p{B{r)) 


— 7r<0i<7r, 


p{9l, - ■ ■ ,9n)'[{ 

i=i 


d9 




( 4 ) 
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where 0 < r < 2y^. For the maximal distance r = 2^/n, the restriction XlILi (^*/^) — 
becomes irrelevant and ii{B{2y/n)) = 1 hy the definition (3). For this reason, we also refer (4) as a 
restricted partition function of the circular unitary ensemhle. Finally, we note that the volume of metric 
halls with other distance measures such as the geodesic distance could he similarly obtained by the 
analytical framework developed in Section III. 

B. Coding-theoretic Bounds 

A fundamental problem in coding theory is to study the maximum size of codewords as a function of 
the minimum distance of a code. A unitary code C with cardinality \C\, 

C = {Ui,U 2 ,...,U|c|} C C/(n) (5) 

is a finite subset of unitary group U{n). With a slight abuse of notation, we define 

r = min |||Uj - UjHf | Ui,Uj € (6) 

as the minimum distance between distinct codewords in U (n). Bounds on codeword cardinality in unitary 
group rely on the normalized volume of the metric ball formulated in (4). In particular, the 

Gilbert-Varshamov lower bound and the Hamming upper bound on the cardinality \C\ are related to 
volumes of metric ball as [8,9] 

h{B{r)) ii{B{r/2)) ' 

'-V-' '-V-' 

Gilbert-Varshamov bound Hamming bound 

The Gilbert-Varshamov bound and the Hamming bound are referred to as a sphere covering bound and a 
sphere packing bound, respectively. The principles of the two coding bounds are different. The Gilbert- 
Varshamov bound (7) is based on a greedy but natural approach to construct a unitary code: start with 
any codeword and keep on adding codewords that have distance at least r from all codewords already 
included until the whole space is covered. Such an algorithm will terminate at \C\ti {B (r)) > 1, where 
{\C\ — l)/i(H(r)) < 1. The argument r/2 in the Hamming bound is related to the error correcting 
capability of the code, where errors made within ‘Hamming sphere’ of radius r/2 can be corrected. This 
generalizes the concept of Hamming distance of q-ary block codes to unitary codes. By taking the union 
of disjoint metric balls each with r/2 packing radius, we arrive at the Hamming bound (7). 

Another important quantity in coding theory, especially in the limit n —)• 00 , is rate of the code [8,9] 

R=-\ogM (8) 

n 
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The bounds on cardinality (7) immediately lead to bounds on code rate 

“ nn ) < ^ < - log2 (— 

Namely, for any n there exists a code in U{n) with a minimum distance r and code rate as bounded 
in (9). The corresponding bounds on minimum distance r as a function of code rate R is obtained as 




( 10 ) 


where denotes inverse function of the volume (4). Note that besides the Frobenius norm other 

metrics such as the spectral norm can be equally considered in the definition of minimum distance (6), 
and the corresponding bounds in (7), (9) and (10) are still valid. 

Finally, we define the diversity sum of a unitary code [2-4], which is closely related to the minimum 
distance (6), as 


S = 



( 11 ) 


The diversity sum (11) is an important performance measure for unitary space-time codes, where a code 
with large diversity sum tends to perform well at the most critical regime of low signal-to-noise ratio [3, 
4]. Therefore, it is interesting to know the largest possible value of diversity sum for a given dimension 
n and cardinality |C|. A useful Hamming-type upper bound to S is given by [4, Eq. (Bl)] 


S < 




( 12 ) 


V n An? 

which is tight for large \C\. Note that the other diversity sum bounds (B2) and (B3) proposed in [4] 
as well as the linear programming bounds [5] can be also analyzed by the framework developed in 
Section III. Finally, from the diversity sum bound (12) and the relation (11), we reach another upper 
bound on minimum distance 

I - 

n4 


r < ^4(^-1 

V n 


(13) 


As discussed in [8], this modified bound is tighter than the minimum distance upper bound (10). 


III. Volume Calculations 

As we have seen, the performance of coding-theoretic bounds (7), (9), (10), (12) and (13) all rely on 
an accurate estimate to volume of the corresponding metric ball. In the following, we derive an exact (29) 
as well as an asymptotic (72) volume formulas using diverse recipes from random matrix theory. Readers 
who are more interested in the applications to coding theory may skip the technical details in this section. 
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A. Exact Volume 

We start by rewriting the n-dimensional integral (4), with the help of a Dirac delta function (5(-), as 


tiiB (r)) = 


^ —7r<0i<7r 


r “ ^ sin^ ^ p 


2 = 1 


dOjdt. 

j=i 


Inserting the Fourier representation of Dirac delta funetion 

<5(t - a) = — / di/ 

J — OO 

into the reformulation (14) and performing the integration over t first, we have 

^ ( 1 - 


1 


(r)) = — 


2tt 


Dnia) do, 


where 


Dn{u) = -/••• 


n 






n 

i=i 


d9j. 


(14) 


(15) 


(16) 


(17) 


Comparing (4) with (16) and (17), we see that the reformulation amounts to eliminating the restriction 
^ r^/4 at the expense of introducing a deformation 11^=1 in the integrand. 

Up to now, the idea of derivation is similar to [4]. To proceed further, we notice that the term 

|2 


n 

l<j<k<n 


^^3 _ 6 ^^^ 


n 

l<_7</!i<n 


qI-Sj _ 


n 

l<j<k<n 


0 _ 0 


= det ) det ) 


(18) 


(19) 


is a product of two Vandermonde determinants, where j,k = 1,... ,n. Invoking Andreief’s identity (see. 
Appendix A), we have 


Dn{o) = - 

C 


r f, 

J det det ^1 (20) 


n\ 


det 




( 21 ) 


Applying the Euler’s formula, the integral inside the above n x n determinant equals the sum of the 
following four integrals 


J cos 6{j — k) cos cos 9^ d9 = 27rcos 

^ J cos0(j —/c) sin cos6*^ d6* = z27rsin^ 

^ J sm9{j — k) cos (^—cos9^ d9 = 0, 

— J sin 9{j — k) sin cos 0^ d9 = 0, 


3 - k 
2 

j - k 


TT 1 Jj—k ( - 


1 Jj—k ( - 


( 22 ) 

(23) 

(24) 

(25) 
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J Y 

^,nj+k+i)j 


n(i) 


X \ 2i+fc 


(26) 


denotes the Bessel function of the first kind [15]. We now have 

Dn{v) = ^ det = det (^Jj_k (0) , (27) 

and hy the definition (26) the Toeplitz determinant (27) has the property 

Dn{- 1 ^) = det {£j) = det (0) = Dn{u). (28) 

Inserting (27) into (16), after some manipulations, we arrive at a one-dimensional integral representation 
for the volume of metric hall (4), 


H{B (r)) 


1 

TT 


sin + sin 



•Jj—k 



du. 


(29) 


The Case n = 2 


In the simplest case n = 2, the integral (29) can he explicitly calculated as 




2,0 

2,2 


1 3 
2 ’ 2 

0,0 


+ G: 


2,1 

3,3 


111 
2 ’ 2 ’ 2 

0 , 0,-2 


(30) 


where a = r^/4 — 1 and 0 < r < 2\/2. Here 


^p,q 


(31) 


(Zl, . . . , CLfi: (^n+1 j ■ ■ ■ j 0,p 
bl, ■ ■ ■ , bfm bjyi-\-lj ■ ■ ■ jbq 

denotes Meijer’s G-function [15,16] (see Appendix B for a brief introduction of this special function). 
The proof of (30) can be found in Appendix C. 


B. Asymptotical Volume 

It seems difficult to explicitly calculate the integral (29) for any n. In addition, from an information 
theoretical point of view, asymptotics of the bounds (7), (9), (10) and (12) are more interesting. These 
facts motivate us to study the asymptotical volume as n goes to infinity. From the representation (16), 
this boils down to the study of the limiting behavior of Dn{a), which is found to be 

lim Dn{v) = e“ 16 . (32) 

n^OQ 

In the following, we will give two distinct proofs of the above strikingly simple limiting formula. The first 
proof is based on matrix-variate hypergeometric functions and the associated zonal polynomials, and the 
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Fig. 1. Dn{u): exact (27) versus asymptotic (32). 


second one utilizes Szego’s strong limit theorem on asymptotics of Toeplitz determinants. Interestingly, 
while the first approach seems to he more instructive and general, it appears to have received far less 
attention in literature. 

Whilst being simple, the convergence of Dn{v) to the asymptotic limit (32) is quite fast as observed 
in Figure 1. As n increases the exact curves oscillate closer to the asymptotic curve and the 

difference is already indistinguishable for n = 4. Note that due to the symmetry (28), we plot only 
positive v in Figure 1. 

Asymptotics of Dnil') via Matrix-variate Hypergeometric Functions 
The equality^ 



allows one to rewrite the n-dimensional integral Dn{v) in (17) as 

Dnia)= [ (34) 

Ju{n) 


*5R{-} denotes the real part of a complex variable. 
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where we utilized the fact that ni<j<A:<n is the Jacohian of the transform (2) from (dU) 

to (ds) n”= dOj. Here, the Haar measure (dU) is normalized to make the total measure unity. The inte¬ 
gral (34) over the unitary group can he represented as a hypergeometric function of matrix argument [17, 
Eq. (91)] as 

[ e*"^{5tdu)}(du) = oFi (35) 

Juin) V 16 y 

The function is referred to as Bessel function of matrix argument [18]. It generalizes the 

classical univariate Bessel function expressed as an integral over the unit circle, cf. the integral inside 
the determinant (21). The general form of matrix-variate hypergeometric function of an n x n Hermitian 
matrix A is represented as [17,19] 


r,Fq (ai,..., Opj 6i,..., bq] A) — ^ ^ ^ ^ 


fc=0 K 


k\ ’ 


(36) 


where the sum over k is a sum over all partitions of integer k into no more than n parts, i.e. k = 


ni + K 2 + ■ ■ ■ + Kn with Ki > k ,2 ■ ■ ■ > Kn > d, and 


(a)« ==n 


{Hj + a- j)\ 


i=i 




(37) 


is the multivariate hypergeometric coefficient [17, Eq. (84)]. In (36), C'k(A) is a zonal polynomial [17, 
19], which is a homogenous symmetric polynomial of degree k in the n eigenvalues of A. Denoting the 
j-th eigenvalue of A hy aj, the zonal polynomial can he represented as [17, Eq. (85)] 


C.(A) = x.(l)x4A), 


(38) 


where 


and 


Xk{1) 


i + j) 


Xk(A) 


det 

det 


(39) 


(40) 


is a Schur polynomial. Schur polynomials form a basis in the space of homogeneous symmetric polyno¬ 
mials in n variables of degree k for all k <n. In particular, it holds [17, Eq. (17)] 

tr'^(A) = J]a(A). (41) 
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In terms of the above notations, we can write DnM as 


Dn{iy) = QFi[n]-—lj 


oo 

(n) 


n i —tcT 

'^K I 16'*'” 


fc=0 K. 


-S)‘ 


^ V ipy ^ 
i.\ 


A:! 

C«(In) 


(42) 

(43) 

(44) 


fc =0 


A:! ^ {n)n 

where the last equality is established by (40). Since the leading order term in {n — j + 1)^. equals 
by the definition (37), for large n we have 


= Win- 3 + 1)«;^. ~ n 
.7 = 1 

Using (41), the sum in (44) for large n becomes 

n) n^oo 


KiH -hKn _ 


E 


in). 


n 


^ J]C,(I„) = 4tr"(In) = l, 


(45) 


(46) 


and we arrive at the claimed result (32), 


OO _ 


lim Dnif) = 

T) —^rv~) f ^ 


16 


k=0 


kl 


= e 16 . 


(47) 


For a detailed account of the asymptotics of matrix variate hypergeometric functions, we refer to [20]. 

We provide here an interpretation of the obtained result. By the identity (33), Dniv) in (34) can be 
understood as the characteristic function of the random variable 




(48) 


i=i 


For any function </(•), random variables of the form Y^^=i9 i^j) called linear (spectral) statistics 
in random matrix theory. Comparing the limiting expression of Dniv) in (32) with the characteristic 
function of a Gaussian random variable J\f ct^). 


(49) 


one observes that the linear statistics (48) follows a Gaussian distribution X J\f {n, fj^) with mean 
/i = 0 and variance = 1/8 as n —)■ oo. Thus, the limiting formula (32) can be interpreted as a central 
limit theorem of the linear statistics (48). Note that this central limit theorem could be also established 
by a group-theoretic approach [21], which did neither use matrix-variate hypergeometric functions nor 
zonal polynomials. 
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Before going to the second proof, we examine the reason behind the fast convergence of Dn{i') to 

1/2 

as n —)• oo. For A = I„, hy repeated use of L’ Hopital’s rule, the Schur function (40) is simplified 


to [22] 


X/ 


/T ^ - i + j) 


As a result, the sum in (44) equals 

C.{ln) 


E 


(n). 


K 

K 


Xk (In) 

{n)K 

i + 3) 


riLi («i + n-j)\ 


For A: < n, hy the orthogonality relation, see e.g. [23, Chap. I. 4] or [21, p. 53], 


(50) 


(51) 

(52) 

(53) 


'^Xk{1)Xk'{ 1) = ^KK'kl 

K, 

we obtain 

E CniJ-nj ^ yj, ^ 


(54) 


(55) 


Inserting the above into (44), we see that the first n terms in the power series of D„(i/) equal the 
corresponding first n terms of - the n —)■ oo limit of Dn{i')- This implies that for any n the 

first 2n + 1 moments of the random variable (48) are exactly the same as the limiting Gaussian random 
variable. Namely, even for n = 2 the first five moments are precisely the same as its Gaussian limit. We 
refer to this fact as finite size mock-Gaussianity. The mock-Gaussian behavior may shed some light on 
the rapid convergence of Dn{y) seen in Figure 1. For n = 2, n = 3, and n = 4, the formal power series 
of Dnii') are computed via (27) as 


i/^ 5i/® 

D 2 M = 1-+- 

^ ^ 16 512 147456 


+ 


113(1/) = 1 - — + E 
^ 16 512 


,6 


24576 


+ 


7z/® 

18874368 

23z/® 

37748736 


+ O , 
+ O ( 1 /^°) , 


(56) 

(57) 


and 


z/^ z/"^ z/® 

DJu) = 1-H- 

^ ^ 16 512 24576 


+ 


1572864 


+ O (z/i°) , 


respectively. Comparing the above with the power series of the limit of 


z/2 Z/^ z/® 

Dooiv) = e 16 =1- 1 - 

^ 16 512 24576 


+ 


1572864 


+ O (z/^°) 


(58) 


(59) 
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one indeed notices that the first 2n + 1 powers of v in the series of match their infinite n 

counterpart (59). Thus, the ‘hulk’ of Dn{y) must he very close to its asymptotic limit with the difference 
mainly resided in the tail as seen in Figure 1. The equality in moments up to such a high order indicates 
that the rate of convergence is quite fast. Indeed, it was proven in [24] that for a class of linear statistics 
including (48), the convergence rate to a Gaussian limit is O for some constant c > 0. This 

super-exponential rate of convergence for linear statistics is uncommon in random matrix theory. 


Asymptotics of Dn{i') via Szego’s Strong Limit Theorem 

As will he seen, the asymptotical formula (32) of Dn{v) is a direct consequence of Szego’s strong 
limit theorem [25] when properly interpreted. We first rewrite the integral (17) as 


Dniiy) = 

/ n 

- 

i2 ^ 

llf{9j)d9,, 

(60) 

C J 



i=i 



where 


f{0) 



(61) 


Following the same steps that led (17) to (21), we have 

Dn{v) = dei{fj-k), j,k = l,...,n, (62) 

where 

= j = 0,±l,±2,..., (63) 

can be interpreted as the Fourier coefficient of the function (61). Thus, Dn{a) is understood as the 
determinant of a Toeplitz matrix {fj-k)'j k=i formed by Fourier coefficients of the function / {9). In this 
form, the limiting behavior of the Toeplitz determinant (62) is known and is given by Szego’s strong 
limit theorem [25], which is stated as follows. If f {9) is integrable over the unit circle and the sum 
E^-oo b'll (ln/)j p converges, where 

(In /), = ^ I" In / {9) d0, j = 0, ±1, ±2,..., (64) 

then 

OO 

\nDn{v) = n(ln/)o + (In/)^ (In/).^ + o(l). (65) 

i=i 

The Landau symbol o(l) is understood as 


lim o(l) = 0, (66) 

n—>oo 


June 25, 2015 


DRAFT 



13 


and in particular (65) is asymptotically tight. In our case there are only two non-zero Fourier coefficients 
for the function In / (9), 

tu 


(In/)- = / cos0e 

47r , 


IV 


= — / cos 0 cos/0 d6* 

dvr 


(67) 

( 68 ) 


(69) 


IV jA j = ±1 

0 / = 0,±2,±3,.... 

Inserting the above into (65) we establish the claimed limiting formula (32). 

Szego’s study of asymptotical Toeplitz determinants [25] was motivated by problems in statistical 
physics in 1940’s. In particular, the spin to spin correlation function of two-dimensional Ising model can 
be written as a Toeplitz determinant [26], whose asymptotics corresponds to the thermodynamics limit of 
such a system. Szego’s original proof [25] imposed strong conditions on f{9), which were successively 
weakened by several mathematicians in subsequent work. The theorem was proven in the present form, 
among others, in [27]. For a concise survey on recent development of asymptotical Toeplitz determinants, 
we refer to [28]. 

We end this subsection by pointing out a by-product of this paper. Since Dn{v) also equals the 
determinant (27), a limiting formula for a Toeplitz determinant of Bessel functions (26) is read off as 

( Jo(x) J-i(x) • • • J-n+l (x) ^ 

Jl (x) Jo (x) • • • J-n +2 (x) 


lim det 

n^oo 


= e 


(70) 


y Jn —1 (^) Jn —2 (^) * * ‘ Jq (^) J 

which seems new. Without using tools from matrix-variate hypergeometric functions or Szego’s strong 
limit theorem, it seems challenging to provide an elementary proof of this interesting asymptotic identity. 


Asymptotic Volume Formula 

We are now in a position to state the main technical result of this paper. Inserting the limiting 
expression (32) into (29), a simple integration yields an asymptotic formula of the volume (4), 


li{B{r)) 


2 /-oo sin ^ -I- sin ( ^ — I ) z/ 
TT 7o 


1 1 / 

= -erf(n)-erf n- 

2 ^ ^ 2 V 2 


-e 16 


dv 


(71) 

(72) 
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Fig. 2. Volumes of metric ball: exact (29) versus asymptotic (72). 


e dt 


(73) 


where 

2 r 

erf(a;) = -^ 

V'^ Jo 

is Gauss error function and (72) is asymptotically tight as n —)■ oo. Due to the fast convergence of Dn{v) 
as observed in Figure 1 and the rapid decreasing nature of the function ^sin ^ + sin 
in (71), we expect that the asymptotic formula (72) approaches the exact one quite fast. Indeed, we 
can see from Figure 2 the approximate volume (72) is already very accurate for n as small as three. In 
Figure 2, the exact volume for n = 2 is calculated hy (30) and the ones for n = 3,4 are obtained by 
numerical integration of (29). 


IV. Applications to Coding Theory 

In this section, we apply the derived volume formulas to the study of coding-theoretic quantities defined 
in Section II-B. As will be shown, the asymptotic volume expression (72) leads to simple yet accurate 
analytical formulas of the bounds^ on cardinality (74), code rate (75) as well as its scaling law (80), 
minimum distance (81), (83) and diversity sum (82). 

^Formally speaking, these are asymptotic bounds valid for n —>■ oo. This fact may not be always mentioned in the following 
as the finite-size accuracy is of primary concern in practice. 
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Inserting the limiting volume formula (72) into (7), the Gilbert-Varshamov lower hound and the 
Hamming upper hound on the cardinality \C\ of a unitary code (5) are obtained as 

2 2 


<|C|< 


erf(n) — erf {n — ' ' erf(n) — erf (n — 

The corresponding bounds on code rate (8) are read off as 


1 


n 


log2 


f(n) — erf (n — y) 


er 


<R<- log2 
n 


f(n) — erf (n — 


er 


(74) 


(75) 


A notably important metric in coding theory is the scaling law of code rate, which establishes existence 
results of codes as dimension grows to infinity [8]. The obtained analytical bounds on code rate (75) 
lead to elementary bounds on limiting code rate. To wit, when dimension n approaches infinity with the 
same speed as the squared minimum distance r^, i.e. under the regime 


n —> oo, r —> oo, with A = — fixed, 

n 


by using an asympfofic expansion of error funcfion 

erf(x) = 1 — 


TTX 


1 + 0 ^ 




X 


oo, 


(76) 


(77) 


and affer some sfraighfforward if fedious manipulafions, if can be shown fhaf fhe following limif exisfs 


lim ^ 

n—^oo Ti^ 


log2 


{\-bf 


erf(n) — erf {n — y) ) (^) 


(78) 


where 


l[a,fe] = 


(79) 


1 A G [a, b] 

0 ofherwise 

defines an indicator funcfion. In (78), fhe limif is faken over fhe regime (76) and In(-) denofes fhe nafural 
logarifhm. As a direcf consequence, fhe normalized limiting code rale R = R/n is bounded by 


(A-2)2 


(A-8)2 




(80) 


41n(2) 

Thus, in fhe asymplolic regime (76) Ihere exisl codes in fhe unilary group wilh code rale per dimension 
bounded by (80). Anolher inlerprelalion of (80) is as a scaling law of code rale bounds, where clearly fhe 
bounds scale linearly wilh dimension in fhe regime (76). Similar resulls for codes in Grassmann manifold 
under a differenl asymplolic regime can be found in [8, Th. 2]. In Figure 3 we plol fhe asymplolic bounds 
on code rale (75) and bounds on fhe corresponding scaling law (80) as a funcfion of minimum dislance 
r. We see lhal as n and r increase fhe rale bounds (75) approach fhe bounds on scaling law (80), as 
expecled. The convergence is reasonably fasl as for n = 16 fhe differences of fhe respeclive bounds 
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Fig. 3. Lower bounds (LB) and upper bounds (UB) on code rate R as a function of minimum distance r: convergence of the 
asymptotic bounds (75) to the hounds on scaling law (80). 


seem indistinguishable for the range of minimum distances r considered even though the scaling law 
bounds (80) are formally valid only when n —)> oo, r —)> oo. It is also observed from Figure 3 that 
the respective lower and upper bounds merge in the high rate regime, whereas their differences increase 
monotonically as the code rate decreases. In fact, from the simple scaling bounds (80) one could easily 
deduce that the maximum gap (at i? = 0) of minimum distances equals for any n. 

By inverting the bounds on code rate (75), the resulting bounds on minimum distance r as a function 
of code rate R is obtained as 

— 2erf“^ (erf (n) — 2^“”-^) < r < 2\j2n — 2erf“^ (erf (n) — 2^“"-^), (81) 

where erf~^(-) is inverse error function. 

We now turn to diversity sum (11), which is relevant to performance of unitary space-time codes. By 
inverting the volume formula (72) and inserting it into (12), after some manipulations we arrive at an 
analytical expression for the upper bound on diversity sum 

S < “ ^erf^ ^erf(n) - . (82) 

The above result makes it possible to analytically study the diversity sum bound. In particular, it can 
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TABLE I 

Diversity Sum Upper Bound (12): Relative Error of the Approximation (82) 


|C| 

24 

48 

64 

80 

100 

120 

128 

1000 

n — 2 

5.0 X 10“^ 

7.3 X 10"^ 

7.8 X 10"^ 

7.9 X 10"® 

7.7 X 10~® 

7.3 X 10"® 

7.1 X 10"® 

1.4 X 10~^ 

n — A 

7.4 X 10“® 

3.4 X lO""* 

5.4 X 10"^ 

6.7 X lO"'^ 

7.6 X 10"^ 

7.9 X lO""* 

7.9 X 10"^ 

2.1 X 10"® 

n = 8 

1.8 X 10“® 

5.4 X 10"® 

9.9 X 10"® 

2.3 X 10"® 

3.5 X 10"® 

4.2 X 10"® 

4.4 X 10"® 

1.2 X lO"’’ 


be verified that the approximative upper bound (82) is a monotonically increasing function of codeword 
dimension n and a monotonically decreasing function of cardinality \C\ with a closed-form relation 
between n and \C\ compactly captured by (82). The exact upper bound on diversity sum (12) was 
numerically evaluated in [4], where the authors stated that “When n is large, the exact computation of 
Vq (the inverse volume (1/|C|)) is rather involved and hence it is also computationally difficult to 
compute the bounds”. Indeed, inverting the integral (29) whose integrand comprises a determinant (27) is 
numerically unstable even for a small n. On the other hand, the fact that the relative approximation error^ 
of the simple analytical formula (82) diminishes quite fast, as illustrated in Table I, makes its usefulness 
more prominent. In Table I, we follow [4, Table I] for different cardinalities \C\ ranging between 24 and 
1000. It is seen that the error of (82) is smaller for not-too-large \C\ than for \C\ = 1000 corresponding 
to the tail of the distribution (72), though the error at which drops rapidly as n increases. Lastly, an 
analytical expression for the modified bound on minimum distance (13) can be similarly obtained by 
inverting the volume formula (72) as 

r < -^\lri^ — (erf“^ (erf(n) — . (83) 

y/n '' 

In Table II we compare minimum distance upper bound in (81) denoted by ri to the modified upper 
bound (83) denoted by r 2 , where the values in each bracket (•, •) are computed from ri and r 2 , respectively. 
As expected, the modified upper bound r 2 is tighter than the upper bound ri. One may also observe that 
the ratio of the two bounds ri/r 2 grows with dimension n, which seems to approach to a constant for 
large n. It turns out that for a fixed code rate R, 

lim — = \/2 lim — = = \/2, (84) 

n^oo r 2 n^oo ^^ (erf(n) - 

^Eor a quantity a and its estimate a, the relative error is defined as absolute value of {a — a)/a. 
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TABLE II 

Comparisons of Minimum Distance Upper Bounds (81) and (83) 



n = 2 

n = 4 

n — 8 

n = 

16 

R = 0.1 

(4.738, 2.588) 

(5.996, 3.969) 

(8.066, 5.656) 

(11.203, 

7.998) 

R = 0.5 

(4.004, 2.828) 

(5.309, 3.971) 

(7.438, 5.605) 

(10.628, 

7.945) 

R = 1 

(3.497, 2.749) 

(4.829, 3.850) 

(6.997, 5.498) 

(10.218, 

7.863) 

R ^5 

(0.802, 0.794) 

(2.251, 2.160) 

(4.912, 4.425) 

(8.379, 

7.138) 

R^IO 

(0.027, 0.027) 

(0.011, 0.011) 

(2.490, 2.429) 

(6.718, 

6.097) 


where the last equality is established hy the asymptotic expansion (77) as well as a corresponding 
expansion of inverse error function 

= T2 \I'" ^ (''' 

Namely, the ratio of hounds tends to a unique constant as n —)■ oo for any code rate R. Though the upper 
hound n is at most \/2 times looser that the modified bound r 2 , their difference is shrinking in the high 
rate regime as seen from Table II. We also see that for a sufficient large n the upper and lower^ bounds on 
minimum distance increase^ with n. This advocates design of large dimensional space-time codes, where 
there potentially exists codes with improved error performance, i.e. with a larger minimum distance. A 
similar conclusion was drawn based on the tabulated minimum distance bounds in [9, p. 3448], where, 
however, some entries of their tables can not be computed due to the small-ball assumption made in [9]. 

V. Conclusion 

In this paper, we derive a limiting volume formula of metric balls in unitary group based on an exact 
integral representation. This simple-to-evaluate asymptotic result is obtained by means of random matrix 
theory. The fast convergence of the derived formula is due to the identified mock-Gaussian behavior. 
The proposed volume formula enables analytical characterization of coding-theoretic bounds of unitary 
codes. Closed-form expressions for bounds on cardinality, code rate including its scaling law, minimum 

''Recall that the lower bound in (81) is always half of the corresponding upper bound. 

^Note that these bounds are not necessarily monotonic increasing functions of n especially in the high rate regime, see e.g. 
the case 7? = 10 in Table II. 
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distance and diversity sum have been derived. These analytical formulas help gain crucial insights into 
the behavior of the bounds. In particular, the largest possible minimum distance gaps between the lower 
and upper bounds as well as between the upper bound and a modified one are quantified. It is also found 
that the code rate bounds obey linear scaling laws in some asymptotic regime of interest. Possible future 
work includes volume computations in other compact classical groups such as orthogonal, symplectic 
groups, and in the Stiefel manifold. 
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Appendix A 

Andreief integral [14] 


For two n X n matrices A(x) and B(x), with the respective ij-th entry being functions Ai{xj) 
and Bi{xj), and a function /(•) such that the integral Ai{x)Bj{x)f{x) dx exists, then the following 
multiple integral can be evaluated as 


Iv 


det (A(x)) det (B(x)) ]^/(xi)dxi = det ( / Ai{x)Bj{x)f{x)dx 


i=l 


( 86 ) 


where V = {a < Xn < ■ ■ ■ < xi < b}. 


Appendix B 
Meijer’s G-function 


The general form of Meijer’s G-function [15,16] is 


^m,n 


Q-l, . . . , 0-12+1? • • • ? Ojp 

^1? • • • ? ^m+l? ' • ‘ tbq 

^ f nr=ir + + ^)n+ir(i-a,-^) 

2vr^ Jc Y{]=n +1 r {aj + z) Y{]=m+i T {I-bj-z) 


X ""dz, 


(87) 


where the contour C is chosen in such a way that the poles of T{bj + z), j = 1,... ,m are separated 
from the poles of T (1 — aj — z), j = 1,... ,n. By definition, Meijer’s G-function is Mellin transform 
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pair of 


I X 


ttl, . . . , Gp 


dx = 


nr=im-+^)n-=ir(i-a,-z) 

nL„+i r(aj + z) ULm +1 r (1 - bj - z) 


( 88 ) 


bi,...,bq I llj=„+l^ -r^;iij=m+l 

In the following, we list some properties and results of Meijer’s G-function that will he utilized in the 
proof of (30): 

• Shifting property 


k^m,n 

kIj p ^ I %Xj 


Inverse argument formula 


/^m,n 

^p,q 


1 


G-l, . . . , Gp 

bl,...,bq 


ttl, . . . , Gp 


= I ^ 


bl,...,bq 

A relation with hypergeometric functions 


= IX 


G\ k^^ Gp “1“ k 
bi + k,... ,bq -\- k 

l-bi,...,l-bq 

1 ( 2 i, . . . , 1 Gp 


pFq (cil ; • • • ? ; • • • ; F) — 


ULiTfe) 

nLir(»0 


Gl'” 


— X 


1 Ul, . . . , 1 (Xp 

0,1 - 6l, ..., 1 - 


An integration identity 


sincx I ax^ 


(Jl, . . . , Gp 

bl,...,bq 


1 _ ^m,n+l 4a 


0, Cli, . . . , dp, 2 

6l, . . . ,6q 


A relation between a Meijer’s G-function and a hypergeometric function 


G^s I ^ 


Ol, 02, O 3 


r (1 - ai + 61) r (1 - ai + 62) 


^ai — 1 


^2; &3 y (r (1 — ai + a2) r (1 — ai + as)) ^ T (ai 63) 

3F2 ^1 — di + bi, 1 — di + 62,1 ~ oi + & 3 ; 1 ~ oi + 02,1 — di + d3; 

A relation between a Meijer’s G-function and a hypergeometric function 


G|° I X 


di, d2 

bi,b2 


= x^^ (1 - (di + d2 - 61 - 62) X 


2F1 (d2 - 62, Ol - b2; di + d2 - 61 - 62; 1 - x) . 


(89) 


(90) 


(91) 


( 92 ) 


(93) 


(94) 
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Appendix C 

Proof of (30): the exact volume for n = 2 
By the definition of Bessel function (26) it holds J_„(x) = {—l)^Jn{x), thus for n = 2, 




the integral (29) can now he split into 


/i(B(r)) = /o --1 +/o(l)+/i --1 +/i(l), 


where 


In{b) = - 

TT 


1 /■°° sinftz/ o /i/ 


- dv. 


Since [15, Eq. 8.442], 

4 (^27 22^r2(n + 1) 1 2 ' 

by the relation to Meijer’s G-function (91), we obtain 


/ 1 

1 F 2 f n + -;n + l,2n + 1;- — 


T 


2 1 

2 


0,0,0 


J = 




8r(3/2) 


^14 A 


0 ,- 1,-2 


By the shifting property (89), the inverse argument formula (90) and the integration identity (92), we 


m = 


-^ / sin6i^ G\'\ — 

271 Jq 4 


1 _i _i 
' 2 ’ 2 ’ 2 


= ( b-^ 


0 , 0,1 


1 _i _i 
' 2 ’ 2 ’ 2 


= —<^33 ^ • 

o,o,-i; 

The case /o(l) is obtained hy using the representation (93) as 

,, , r 2 (i/ 2 ) ^ /I 1 \ 1 

-^o(l) - —— 31^2 ( t:, t:, 0 ; 1 , 1 ; 1 ) - 
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The integral Ii{b) is calculated by following the procedure of Io{h) as 


1 /’°° / 

= wifi, Ou t 


0 

1 _i _3 

2 ’ 2 ’ 2 




—( 6-2 
27r 3,3 


0,0, i 

i _1 _3 
2 ’ 2 ’ 2 


27r^3,3 


= —G: 

27r 


2,0 

2,2 


_i i 3 

2 ’ 2 ’ 2 

0 , 0 , 

1 3 
2 ’ 2 

0,0 


(105) 


(106) 


(107) 


(108) 


where the last equality is established by the definition of Meijer’s G-function (87). The case /i(l) is 
obtained by using the representation (94) as 


7i(l) = 0. (109) 

Inserting (103), (104), (108) and (109) into (96), we obtain the claimed result (30). 
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